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The parametrically driven Ginsburg-Landau equation has well-known stationary solutions — the 
so-called Bloch and Neel, or Ising, walls. In this paper, we construct an explicit stationary solution 
describing a bound state of two walls. We also demonstrate that stationary complexes of more than 
two walls do not exist. 
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I. INTRODUCTION 

In this paper we derive a new class of exphcit and 
physically meaningful solutions to an equation that has 
been under scrutiny, in various contexts, for more than 40 
years. The equation is the resonantly driven Ginsburg- 
Landau; in its most general form it reads: 

= ifi+iiy)iJ+igi+ici)i!^^-{g3+ic3)\^\^^-h{^*)''-\ 

(1) 

Eq.(^) describes a one-dimensional chain of coupled 
weakly nonlinear self-sustained oscillators in the contin- 
uum approximation. The chain is subjected to the peri- 
odic forcing at the frequency « ncuo, where ujq is the 
frequency of the undriven spatially homogeneous linear 
oscillations. The complex variable ip = ^{x, t) is the 
slowly varying amplitude of the resulting nonlinear oscil- 
lations with the period T = 27rn/r2. 

Before specialising Eq.(^) to the particular case that 
will concern us in this paper, we briefly comment on the 
physical meaning of its coefficients. First of all, the pa- 
rameter /i > measures the distance to the supercriti- 
cal Hopf bifurcation at which a (spatially homogeneous) 
stable limit cycle appears. For the stability of the self- 
sustained oscillation one also needs 53 > 0. The real 
constant v is proportional to the frequency detuning, 
V (X Q,/n ~ LjQ. In the derivation of J^l j it is assumed 
that both the linear growth rate /i and the detuning v 
are small compared to the frequency ujq. (See e.g. [l|.) 
Next, the parameter C3 describes the nonlinear frequency 
shift, while the derivative term accounts for the interac- 
tion of the neighbouring oscillators in the chain, with the 
real and imaginary part of the coefficient gi + ici per- 
taining to the dissipative and reactive type of coupling, 
respectively. Finally, h is proportional to the amplitude 
of the forcing. The (n — l)-st power of -0* in Eq.|p[J re- 
sults from the resonance forcing of order n : 1 Ql- The 
1 : 1, 2 : 1, 3 : 1 and 4 : 1 resonances have been studied 
most extensively in literature. In this work, we focus on 
the case n = 2. 



The analyses of Eq.(^ usually start with considering 
the gradient, or variational, limit 0,01 ■ The variational 
limit corresponds to the assumption that the coupling 
is purely dissipative (ci = 0), frequency detuning v is 
zero and the homogeneous oscillations are isochronous 
(c3 — 0). With these assumptions, a suitable rescaling of 
t, x, V' and h produces 



IV'lv + ^/'-^v'* 



(2) 



In this paper, we consider stationary solutions of 
Eq.(|2Jl; these satisfy 



-'4'xx-\i^?'>P + i'-H* 



0. 



(3) 
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We will take h to be positive in what follows; this can 
always be achieved by an appropriate phase-shift of -0. 
Historically, Eq.© was first introduced in the context of 
the anisotropic XY-vaodel, which was used to model an 
easy-axis ferromagnet near the Curie temperature ,5, ^. 
Nonstationary magnetization configurations were consid- 
ered as solutions to Eq. Q f?! . The investigations of the 
more general Eq.QJ (still with n — 2, though), includ- 
inganalyses of small nonvariational effects, were reported 
in |1, 11,01 ■ Apart from the magnetic applications, these 
studies were motivated by research in liquid crystals [10| , 
experiments with the periodically forced light-sensitive 
Belousov-Zhabotinsky reaction [Tl| and work done in op- 
tics, in particular with regard to optical parametric oscil- 
lators '12l| and lasers with intracavity parametric amplifi- 
cation JJJ. Eq.QJ also appeared as a phenomenological 
equation for t he p arametrically excited surface waves in 
viscous fluids |14| and granular media Iisj . 

The nontrivial (spatially nonhomogeneous) solutions 
which attracted interest in the context of each of these 
fields are domain walls, or kinks, also known as dark 
solitons in nonlinear optics. The domain wall is a 
localised interface between two different homogeneous 
backgrounds (known as domains in the magnetic con- 
text). The possible backgrounds are described by the 
constant nonzero solutions of Eq.©, -0 = and 

= ±iA+, where A± = ^/l±h. The first of these is 
known to be unstable while the second is stable. We will 
only consider dark solitons propagating over the stable 
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known explicitly, for /i = [lMl20l|: 



_(c) 




FIG. 1: The Neel wall (a), and the left- and right-handed Bloch 
walls [(b) and (c) respectively]. The solid line corresponds to 
the real part, and the dotted line to the imaginary part. Here 
h = 0.05; we plot the two walls for a small value of the parameter 
in order to accentuate the difference in their widths. 



background — asymptotically, all of the solutions con- 
sidered here will satisfy IV'P ^ ^+ sls \x\ oo. 

The soliton solutions to Eq. with the desired asymp- 
totic behaviour can either be topological (with a phase 
difference of 180° between the two asymptotic values) 
or nontopological (with no change in phase between the 
asymptotic fields). Eq.|(2Jl admits two explicit topological 
solutions. The first is the Neel wall [sL ll ITl[l7| . 



ip^ix) — ij4+ tanh(^_|_a;). 



(4) 



so named because its magnitude \ip\ vanishes at the ori- 
gin, at which point the phase becomes discontinuous. (In 
magnetism, a Neel wall is a domain wall with vanishing 
magnitude of the magnetisation vector at its centre.) ft 
is also called the Ising wall, the name appealing to the 
Ising model which was used to emulate the easy-axis fer- 
romagnet. The Neel wall exists for all positive h. 
The second topological solution has the form 



iPb{x) — iA+ tanh(i?a;) ± Csech(i?a;) 



(5) 



where B = and C = ^1 - 3h |E 111 El- This 

solution is usually referred to as a Bloch wall (which, 
in magnetism, is a domain wall connecting the two do- 
mains smoothly, with the magnetisation vector remaining 
nonzero everywhere). The Bloch wall exists in two chi- 
ralities, distinguished by the sign of the real part in ((SJ. 
For convenience, we will refer to the solution with posi- 
tive (negative) real part as the left-handed (right-handed) 
Bloch wall. Regardless of the chirality, the Bloch walls 
only exist for h < ^. 

In addition to the Bloch and Neel walls, Eq.© pos- 
sesses nontopological solitons. One such solution is 



3 n 

1 - -sedi{Bx) 



SBtaiih{Bx)sech{Bx). (6) 



Here B = y/ih = ^^4/15. In addition, a class of bubble- 
like solutions (including Eq. © as a particular case) was 
found numerically [23. As is usual with numerical so- 
lutions, the structure of the solitonic bubbles has not 
been completely understood. It has also remained un- 
clear whether they form a continuous family (families) 
and if they do, what is their range of existence (in h). 
Finally, a pertinent question is of the generality of these 
solutions: can the resonantly driven Ginsburg-Landau 
equation support localised, stationary structures other 
than the walls and bubbles? The aim of this paper is to 
address all of these issues. 

Here, for each /i < i, we construct a two-parameter 
family of explicit bubble-like solutions and demonstrate 
that they can be interpreted as stationary bound states, 
or complexes, of a Bloch and Neel wall. Physically, the 
constructed solutions represent phase domains of arbi- 
trary length, which are 180° out of phase with the back- 
ground field. We also prove that these solutions exhaust 
the list of possible stationary complexes — there can be 
no stationary bound states of more than two walls. 

It is appropriate to mention here that equation Q ap- 
peared previously in several other, unrelated, physical 
contexts. Written in terms of the real and imaginary 
parts of it coincides with the stationar y equat ion for 
the so-called MSTB model of field theory [Illll. More 
recently, it occurred as a stationary limit of the paramet- 
rically driven nonlinear Schrodinger |2l| equation. The 
latter equation arises in a large variety of physical appli- 
cations including nearly resonant optical parametric os- 
cillators |22| . easy-plane ferromagnets in magnetic fields 
[2ll |. and surface waves in wide, vertically vibrated chan- 
nels of inviscid fluid 0, . Accordingly, we expect the 
new solutions to admit physical interpretations in these 
fields as well. 

The paper is organised as follows. In sectionniwe de- 
rive a family of explicit two-soliton solutions and cast it 
in a symmetric form allowing easy visualisation. These 
solutions will be interpreted as bound states of a Bloch 
and Neel wall. Several families of singular solutions ap- 
pear as by-products in this construction; these will be 
used later for auxiliary purposes. After that (section 
IIII|) we show that there exist no other bounded solutions 
of Eq. ^ which would asymptotically approach the sta- 
ble {zLiA+) background. This means that apart from 
the constant solution, the only nonsingular solutions to 
Eq.© are the Bloch and Neel walls, and the Bloch-Neel 
complexes constructed in this paper. Finally, some con- 
cluding remarks are made in section Hvl 
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II. EXACT SOLUTIONS FOR THE 
BLOCH-NEEL BOUND STATE 

A. The Hirota construction 

In order to construct a bound state of two walls ex- 
plicitly, we employ the Hirota bilinear formalism (see e.g. 
111). Letting 



where G is complex and F real function of equation 
is cast in the bilinear form: 

F[DlG • F + (2 - \)GF - 2hG*F] 
-G[DlF ■ F - XF^ + 2\G\^] = 0. (8) 

Here Dx is the Hirota operator defined on ordered prod- 
ucts of functions: 

D^a-b={dx-dyra{x)b{y)\x=y. 

In ((SJ we have added the term XGF^ to the second 
bracket and subtracted it from the first one. The con- 
stant A will be chosen later. 

By making the substitution ((7|) we increased the num- 
ber of unknowns while the number of equations remained 
unchanged. We can use this freedom to set the first and 
the second term in © to zero (separately): 

Dlu-F+(2Al~X)uF^0, (9) 



We now choose A = 2^4^. Then, equations (fT^ - lfl^ give 
vq = A. (Here, and in the rest of the paper, A stands for 
A^.) Next, at the order e-^ we obtain 

i-dl + Ah)ui = 0, (14) 



dUAFi+vi)=0, (15) 



{~dl + 2A^)Fi^2Avi. (16) 

Equation yields ui ~ e^^, where Q\ = 2h^/'^{x — xi) 
and xi is an arbitrary constant. From H15|l we infer vi — 
—AFi. Substituting this into Hlt)|) we get 

{-dl + AA^)Fi^Q, 

whence Fi = e^^, where 62 — 2A{x—X2) and X2 is another 
arbitrary constant. 

The order gives three equations: 

{-Dl + Ah){u2-l + ui-Fi)^Q, (17) 



Dl{vQ-F2 + vi-F^+V2-l)^0, (18) 

{-Dl + 2A^){Fi- Fi + 2F2-l) = 2{ul + vl + 2vQV2). (19) 
Substituting for ui and Fi in 1171) this equation becomes 



DlvF+ {2A\ - X)vF = 0, 



DlF-F-XF^ + 2{u^ + v^) = 0, 



(10) 



(11) 



where we let G = u + w and decomposed the first bracket 
in (jSl into its real and imaginary part. Now we look for 
a solution to the system j^- pi|) as a series 



U = £1*1 + e U2 + 



= fo + ei'i + e 1^2 + 



F = 1 + eFi + £2^2 + • • • , 

where e is a formal expansion parameter. An explicit 
solution will arise if the series truncates at a finite power 
of e. 

Substituting into ll^- llll) . the order ^ gives 



a2«o + (2^^-A)i;o = 0, 



(12) 



{^dl + 4/i)u2 = 4AC_e''i+''^ , (20) 

where 

C-=A~2h^l'^. (21) 
Ignoring its homogeneous solution, we get 



A - 2\/7i g 

U2 = !=& 

A + 2Vh 



J1+O2 



C4 



(22) 



Next, substituting vq = A and vi — ~AFi into H18|l . we 
obtain V2 — —AF2, after which equation 119|l becomes 

{-dl + 4A^)F2^e^'\ 

whence, ignoring again the homogeneous solution, 

1 



F2 = 



4C+C- 



(23) 



The singularity arising for A^ — Ah can be removed by 
letting xi = 00 in 01. 

To the cubic order in e we get 



2v?:~X^ 0. 



(13) 



i-Dl + Ah){ui ■ F2 + U2 ■ Fi + U3 ■ 1) ^ 0, (24) 
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Dlivo ■F3 + vi-F2+V2-Fi+V3-l) = 0, (25) 



(-D2+2A2)(F3-1+F2-Fi) 

= 2(uiU2 + t;iW2 + V0V3). 



(26) 



Substituting for wi, M2, -Fi and i^2, equation 1)24(1 becomes 

i-dl + Ah)u3 = 0, 

whence U3 = 0. On the other hand, equations 
reduce to the system 



{~dl + AA^)iv3-AFs)=0, 
which gives U3 = AF3 and 

C- 



F. 



^201+92 



Next, the order yields 

{-Dl + 4/l)(u4 ■1+U2- F2+U1- F3) ^0, 



(27) 



(28) 



(29) 



(30) 



Dlivi ■l + V3-Fi+V2-F2+vi-F3+vo- F4) = 0, (31) 



i~Dl + 2A^){2Fi • 1 + 2F3 • i^i + F^) = 

2{ul + V2+2viV3 + 2voV4), (32) 

where we have taken into account that U3 — 0. Substi- 
tuting for all the variables in (|30|l we obtain 

{~dl + 4h)u4 = 0, 

whence U4 = 0. Equations l|31ll - 132() reduce to the homo- 
geneous system 

dl{v4 + AF4)^0, 
{-dl -f 2A^)F4 = 2Av4. 

We choose the trivial solution, V4 — F4 — 0. 

The order is the last order that we have to do "by 
hand" ; in dealing with all higher orders (e" with n > 
6) we will simply invoke the machinery of mathematical 
induction. To e^, we get 

i-Dl + Ah){u2-F3+U5-l)=0, (33) 



Dlivo ■F5+V2-F3+V3-F2+V5-l) = 0, (34) 



i~Dl + 2A^){F5 • 1 + F2 • fg) = HAv^ + V2V3), (35) 

where we have taken into account that U3 — U4 ^ V4 ~ 
F4 = 0. Substituting for U2 and ^3, equation H33|l gives 

i-dl + 4h)u5 = 0, 

whence — 0. Recalling that V2 = —AF2 and — AF3, 
equation becomes 

dUv5+AF5)^0, 

whereas making use of {—D^ + 2A^)F2 ■ F3 — 2V2V3 in 
(|^ we get 

{^dl + 2A^)F5 - 2Av5 ^ 0. 

The last two equations are satisfied by letting V5 = -F5 = 
0. 

Finally, we prove that all coefficients u„,u„, and Fn 
with n > 6 are also equal to zero. We assume that u„_i — 
Vn-i — F„-i = and show that this entails u„ = v„ — 
Fn — 0. Consider, first, equations ((Tn)l -l(TT |l . Setting to 
zero the coefficients of e", we obtain 



Di 



VQ ■ Fn + y^^Vk ■ Fn 



0, 



(36) 



k=l 



{-Dl + 2A^) Fk ■ Fn-k + 2Fn ■ 1 j 

n-l 

= 2'^{ukUn-k + VkVn-k) + 4voVn- (37) 
k=l 

Since Un-k = for all3<rt — fc<ri— 1 and Fn-k — 
Vn-k — for &\\ A < n — k < n — 1^ the sum involving 
Un-k iu the right-hand side of H37|) begins with k = n—2 
(rather than with k — 1), while all sums involving Fn-k 
and Vn-k begin with k = n — i. On the other hand, 
since n — 2 > 4, all Uk in the sum in (|37|l arc equal to 
zero. In a similar way, all Vk and Fk in the sums in l|36|) 
and (|37|l equal zero - except Vn-3 and F„_3 for ti = 6. 
Therefore, for n > 7 equations (|36|) - (|37|l become a pair 
of homogeneous equations for Vn and Fn] hence we can 
set Vn = -Fn = 0. For n = 6, we get 

Dl{vo-Fe+V3-F3 + ve-l)=0, 
{-Dl + 2A^){F3 ■ F3 + 2Fe ■ 1) = 2{vl + 2voVe). 

Using V3 = ^^3, this also reduces to a homogeneous 
system for vq and Fq. 

Lastly, equation ^ gives, to the order e": 



{-Dl + Ah) 



' Uk ■ Fn-k +Un-l \ = 0. 



Since n — 3 > 3, we have Uk — for fc > n — 3. On the 
other hand, all Fn-k = for A: < n — 3 and so all terms 
in the sum equal zero. Hence u„ = 0. Q.E.D. 



5 



B. The explicit solution and its interpretation 

Thus we have constructed an expUcit solution of the 
form i/j — {u + iv)F~^ , where u,v and F are polynomials 
of e^^ and e^^ with real coefficients. It will be shown 
later that if C_ < 0, the solution has a singularity. [We 
remind that C_ is given by Ea. (|21|l .] For now, we assume 
that C_ > and cast the solution in a more symmetric 
form. 

First of all, the parameter e can be absorbed into e^^ 
and e^^ through the re-definition of the arbitrary con- 
stants xi and X2- Next, we define xi = 9i — a — P and 
X2 = ^2 - 2/3, where 



= 4C+C_ = 4(^^ 



A 



-4/i), 



(38) 



The new phases xi and X2 still involve arbitrary con- 
stants xi and X2. We can choose these constants in such 
a way that, up to an overall translation, 



Xi 



2/ii/2(x-s), X2^2A{x + s), 



(39a) 



where s is the only free parameter remaining. The nu- 
merator and denominator of the solution 



Tp — {u + iv)F ^ 
are then given by the following expressions: 



V = A{1 - e 



2/3+X2 



e 

o2/3+2xi 



o2xi-(-X2^ 



(39b) 



(39c) 
(39d) 
(39e) 



Asymptotically, ^ lA both as a; — > oo and x — oo, 
and hence Ea. H39|) has the form of a bubble. If we per- 
form the reflection x —^ —x, and at the same time replace 
s with — s, the real and imaginary parts of the solution 
change according to u/F ^ —u/F, v/F — > v/F. Denot- 
ing the solution (|39() by ?/'(x;s), we therefore have the 
following symmetry: 



^(- 



-x: 



'ip*{x;s) 



This implies, in particular, that the solution H39|l with 
s = has an odd real and even imaginary part. When 
h — j^, the s = solution reproduces the explicit solu- 
tion that has been known before, Eq.®. (More precisely, 
it is equivalent to © with x — > —x.) 

The solution 1)39(1 describes a stationary complex of a 
Bloch and a Neel wall, with the parameter s character- 
ising their separation. This is easily seen by examining 
Ea. (|39|l in the limit of large s. First, let s be large and 
positive. For a; ~ s, we have e^^ ^ e^^^ ^1. In this 
region, the solution reduces to 

ip = lAtanhX — Vl — 3/i sechX, 



(a) 



.■ h = 0.2 
■■ s = -10 





} 


ir- 


(b) 


h = 0.2 

s = 













(c) 




r 

h = 0.2 





(d) 



h = 0.2 
s = 10 

— 1^* 



FIG. 2: The bubble solution ^ for (a) s = -10, (b) s = and 
(c) s = 10. The panel (d) corresponds to the bubble of opposite 
chirality [i.e. with tp{x) —tp*{x)]; like (c), it is plotted for 
s — 10. The solid line corresponds to the real part, and the 
dotted line to the imaginary part. 



with X = 2\fh[x — s) — (3, which is a right-handed Bloch 
wall centered at xq — s + -^ph'^^^ . In the region x ~ — s, 
we find that (|39|l becomes 

i(; = -iAta.nh[A{x + s) + /]], 

which is a Neel wall centered at xq = —s — PA~-^ . On the 
other hand, if we let s be large and negative, then we find 
a Neel wall on the right (centered at xq = —s + f3A~^) 
and a right-handed Bloch wall on the left (centered at 
xo = s - i/3/i-i/2). 

These conclusions are illustrated by FigI21 which de- 
picts the real and imaginary part of H39|l for several rep- 
resentative values of s. We also note that the transforma- 
tion X —^ —X, s ~* ~s which was noted above to change 
the sign of the real part while leaving the imaginary part 
intact, simply flips the chirality of the Bloch wall. (See 
FigHd)). 

There is yet another way of seeing that Eq. (|39|l repre- 
sents a bound state of two walls, and this time we do not 
have to assume that s is large. Consider the integral 



{A^ - IV^H dx 



(40a) 



which gives an integral measure ( "area" ) of solutions with 
|?/'(x)| — > A as |a;| ^ oo. This integral usually has some 
physical meaning in the nonlinear Schrodinger-based in- 
terpretations of Eq. (PI) . For example, if Eq. ^ is regarded 
as a stationary reduction of the Landau-Lifshitz equation 
for the easy-plane ferromagnet in the external magnetic 
field llj, the integral (|lla|l gives the total number of 
magnons in the excited state of the ferromagnet. In the 
Ginsburg-Landau based applications of Eq.©, the in- 
tegral (|40a|l does not usually have any special meaning. 
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but nevertheless can be used as a scalar characteristic of 
stationary solutions. Letting ij; = {u + iv)F~^, Ea. (|40a|) 
becomes 

1 = I (A^F'-u'-v') ^. (40b) 

For the Bloch wall, this integral equals Ib = 2B while 
for the Neel wall we get /n = 2A. As for the solutions 
(|39|l . the evaluation of the integral (|40|) becomes trivial 
if we notice that, due to Ea. (|ll|l . the integrand is a total 
derivative: 

{A^F^ - u2 - 1,2) ^ (d^F/F) . (41) 

Evaluating dxF/F at a; = ±oo we conclude that each 
member of the family of bubbles H39I) has the "area" / = 
2{A + B) — which is exactly the sum of Jb and Jn- 

We conclude this section by mentioning that bound 
states of a pair of stationary dark solitons have been also 
found in the (integrable) Manakov system . 

C. Other solutions 



If the top sign is chosen in H44() . an argument similar to 
the one following Eq. (|42l) demonstrates that F must pass 
through zero whereas u is strictly positive; once again, 
u/F is singular. If the bottom sign is chosen, F is positive 
in the asymptotic regions x — *■ ±cxd, and this type of 
argument would not work. However, at the point x — 
— s — /3A~^, we have X2 = — 2/3 and so Eg . (|44c|) (with the 
bottom sign chosen) becomes 

F{x) = e^^' {e-^^ - e^^) . (45) 

By Ea. H43|) . e^^ > 1, which means that the contents of 
the brackets in H45|) is negative. Thus for this value of 
X, F < 0. But since F is positive in both asymptotic 
regions and is continuous, it follows that F must vanish 
at least at two points. Since the bottom sign has been 
chosen, (j44bp implies that v > A > 0, so when F = 0, 
the quotient v/F is infinite. 

Thus, in summary, of all the families of solutions found 
in this section, only family H39|) is regular, and only for 
/i < i. The other solutions are singular and not physi- 
cally relevant. 



If we had chosen Fi = — instead of at the order 
e, we would have arrived at a different solution. Assum- 
ing /i < i (so that C_ > 0) and defining a, /3, xi and X2 



as m Eas. l|5S|l and H39a|) 
as Ea. (|39bjl with 



this new solution can be written 



A{1 + e2'^+^^ 



e 



(1 



F = 1 



g2/3+2xi _ g2xi+X2^ 



(42a) 
(42b) 
(42c) 



It is not difficult to check that this solution is singular. 
Indeed, the function F{x) is continuous, and has opposite 
signs at the two infinities: -F ^ 1 as a; ^ — oo and F — s- 
— oo as a; oo. Therefore, it must pass through zero at 
least once. Since u cannot vanish for finite x, we have 
u/F = 00 wherever F = 0. Although singular solutions 
are not physically meaningful, there will be some indirect 
use for them in the next section. 

Assume now h >h- Here we have C_ = A — 2h^^^ < 
and the functions (|39|l and (|42|l are no longer solutions, 
because the parameters a and /?, as defined by Ea. (|38|l . 
are not real. We can, however, define a and (3 by 



= -4C+C_ =4(4/1-^2), 



C- 



2/ti/2 ^ A 
2/1I/2 



(43) 



and, instead of Eas. l|39c(l - I|39e|l and l|42a|l - H42c|l . arrive at 
the following two solutions: 



u 

v = A{\T 6^'^+^=^ 4 
F = l± 6^1"+^' 



e 

o2/3+2x 



(l±e>^=^), 

2X1+X2N 



g2/3+2xi -p g2xi+X2 



(44a) 
(44b) 
(44c) 



III. "COMPLETENESS" OF THE LIST OF 
SOLUTIONS 

It is convenient to rescale Eq. (PJ so that the asymptotic 
values are equal to ±1. Letting ipix) ~ iA'i>{Ax), Eq.Q 
becomes 



A2 



21^-12^ 



2 1 



A2 I * 



(46) 



We will demonstrate that solutions ^ and l|39|l 
(when suitably rescaled) are the only bounded solutions 
of Eq. (|46|1 that asymptotically approach one of the stable 
fiat backgrounds \E' = ±lasa:^— 00. We will consider 
solutions of Eg. 1)46(1 as trajectories in a 4-dimensional 
phase space and show that for any direction in which 
a trajectory can leave the fixed point (representing the 
fiat background), we already have a solution leaving in 
that direction. By uniqueness, no other trajectories will 
be allowed to exist. For definiteness, we confine our- 
selves to the case of trajectories that approach ^ = 1 
as X — !■ — 00; a similar argument is valid for solutions 
approaching ^' = — 1. 

For convenience of presentation, we start by listing all 
known solutions to the rescaled equation H46|l . Besides 
the solutions derived in the previous section, we also in- 
clude solutions which are the unbounded counterparts of 
the isolated Bloch and Neel walls. In all of these, we 
make explicit the translational invariance of Eg. (|46|l . 

For /i < |, apart from the flat background ^Pq = 1, the 
list consists of the Neel wall = — tanh(a; — .tq) and 
its singular counterpart 



' coth(x — xq); 



(47) 
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the Bloch walls 

= - tanh [B{x - xq)] ± iCsech [B{x - xq)] , (48) 

where B = 2{A^ - l^/^A'^ and C = (4 - 3^2)1/2^-1^ 
and, finally, the Bloch-Neel complex (the bubble) which 
we write together with its unbounded counterpart: 



u + iav 



BN 



I ± g2/3+x2 _|_ e2/3+2xi ± e2xi+X2 ■ 



(49a) 



Here 



=2/3 _ 



= 1 ^ e^'^+x^ - e^^+^xi ± e^x^+x^ , (49b) 
V ^2{1 + B)e^^lTe^^), (49c) 

= (1 + B){1 - B)-\ XI = B{x ~ xo ~ s) and 
X2 = 2(a; — Xq + s). The real parameters s and a;o are 
arbitrary; s characterises the separation distance and xq 
describes uniform translations. The sign factor cr = ±1 
determines the chirality of the Bloch wall bound in the 
complex. [Note that <j is not correlated with the sign fac- 
tor in (|49b|l - H49c|l .] If the top sign is chosen, Eq.lg^ gives 
the rescaled regular solution (|39|) . If the bottom sign is 
chosen, we have the rescaled singular solution (|42|l . 

For /i > i, the functions 'So, ^'n and defined in 
the previous paragraph remain as solutions. Eas. (|48|l and 
(|^ . on the other hand, are no longer solutions. Eq.l^HJ 
is replaced by 

= - coth [B{x - xa)] ± zCcosech [B{x - xq)] , (50) 

where B stays as it was previously defined whereas C is 
now given by C = (SA^ - 4)^/2^-1 . Eq.lgHl is replaced 

by 



4' 



u + iav 



BN 



1 zp e2/3+X2 _ e2/3+2xi ± e2Xl+X2 ' 



(51a) 



where 



w = 2(1 + 5)6^^1(1 Te^^), (51c) 

and all parameters are defined as for (|49|l except e^^^ — 
{B + 1){B - . Both (inni and l^ll) are singular solu- 
tions [for all choices of signs in (|^ ] . 

We now turn to Ea. H46|) and write it as a dynamical 
system for a particle on the plane: 

r/xx-2(r/2+e2)r/ + 2r/ = 0, (52a) 
e..-2(r?2+ ^2)^^(2-52)^^0, (52b) 

where 5* = 77 + with 77 and ^ real, and B^ = 4(^2 „ 
l)v4^2 (as defined for Eq.l^SJ). This is a hamiltonian 
system, with the hamiltonian 

n^lw.+ii - +e - If - B^e]. (53) 

There is also a second, independent, integral of motion: 

I = (Cry. - v^.? + H - v'e - iv' - m . (54) 



(In the derivation of (|54() we were guided by the results 
of [23 which considers a similar system.) Thus all tra- 
jectories of the system (j52|l are confined to lie on a 2- 
dimensional surface, defined by the constraints (|53|) and 

(EH). 

We consider trajectories that flow out of the fixed point 
iV: Vx,^x) = (1, 0, 0, 0). (On these trajectories, the con- 
served quantities obey Ti. = = X). This fixed point is a 
saddle, with two positive and two negative real eigenval- 
ues. Thus in a neighbourhood of the fixed point on the 
unstable manifold, rjj. will be either positive or negative 
(as will ^2,), while oscillatory behaviour is not possible. 
In this neighbourhood, the sign of rj^ will obviously be 
the same as the sign of {rj — 1); similarly, the sign of 
will be the same as that of (^ — 0). As for the magnitudes 
of rjx and these are determined by rj and ^ via the 
constraints H — 0, T — 0. Thus the variables 77 and ^ 
uniquely determine rjx and in the vicinity of the fixed 
point and therefore provide coordinates on the local un- 
stable manifold. 

In the vicinity of the fixed point, Eas. H52a|l and (|52bl) 
can be written as 



4(77-1) + 2^2^ 
Cx. = B^C 



(55a) 
(55b) 



(Note that we cannot drop the ^2 term in the top equa- 
tion here, as the equation does not include any term lin- 
ear in ^.) The solution to (|55|l satisfying 7; — > 1, ^ — > 
as a; ^ —00, is 



77 = 1 + Xe2 



AA2 



2Bx 



2(1-^2) 



(56a) 



(56b) 



The real constants M and Af are arbitrary and provide a 
parametrisation of the local unstable manifold: each pair 
{Ai, J\f) defines a trajectory on the manifold, and the 
other way around — for each point (77, ^, 77^;, ^x) on the 
manifold, we can find a pair (A^, Af) such that there is a 
trajectory connecting (77, ^, r/x, ^x) to (1, 0, 0, 0) which 
is described by H56|l . Indeed, for the given pair of coor- 
dinates (770, Co) we can solve to get 



7We2"° = 770 



eo 



2(1-^2) 



-1, 



j^^Bxo 



Co. (57) 



Due to the translational invariance of the system H52|l . 
replacing x with x — xq in (|56|l simply furnishes a dif- 
ferent parametrisation of the same trajectory. Hence the 
required pair (A^, Af) is obtained, e.g. by setting xo = 
in |57|). 

Thus, if we factor the translation invariance out, there 
is a one-to-one correspondence between {M, Af) and 
(r/, ^) on the local unstable manifold. We now show 
that for every pair {A4, Af) (with —00 < Ai, Af < 00) 
we already have an explicit solution in our list, regular 
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or singular, with the asymptotics (|56|l . This will imply 
that our list is complete and no other solutions with this 
asymptotic behaviour can exist. 

First of all, the solution corresponding to {A4, A/") = 
(0, 0) is the flat background = 1- Keeping Af — 
0, we have two possibilities (for all h): in the case 
Ai > 0, Eas. (|56|) give the asymptotics of the solution 
"I'm, Ea. l|T7|l . with xq defined by = 26"^"^°; similarly, 
the case M < corresponds to the Neel wall ^^n, with 
M = -2e'2^o. 

Now let = while Af ^ 0. Here the two cases, 
/i < i and h > ^, have to be considered separately. For 
/i < i, the pair (0, Af) corresponds to the Bloch wall 
(EHI), with 

|AA| = 2(1 - 3hy^^A-^e-^'">. (58) 

The sign of Af is arbitrary and determines the chirality 
of the soliton. If ft, > |, we recover the solution ^b', this 
solution also occurs with two different chiralities, depend- 
ing on the sign of A/". 

Finally, we let A4 ^ 0, A/" 7^ 0. Assume, first, that 
ft < i and consider the solution Ea. l|l^ . This 

solution comes in two chiralities, a — ±1. We let Af > 0; 
this selects one of the chiralities {a = +1). (The case 
Af < can be considered in a similar way.) Comparing 
the asymptotic behaviour of the solution (|49|l to Eas. H56|l . 
we find 

A4 = zp2e^('^""=°+'') 

and 

Af =2{1 + B)e-^^''"+'\ 

Thus A4 < corresponds to the (regular) Bloch-Neel 
bound state [top sign in 1311)], and A/{ > to its singular 
counterpart [bottom sign in 149|l ]. For ft > i, a similar 
consideration involving Ea. (|51|l demonstrates that any 
pair {A4, A/") with A4 ^ 0, JV ^ corresponds to a known 
solution as well. Q.E.D. 

IV. CONCLUDING REMARKS 

In this paper, we have explicitly constructed bound 
states of a pair of domain walls in the resonantly forced 
Ginsburg-Landau equation. The constructed solutions 
represent domains (of arbitrary length) that are 180° 
out of phase with the background field. We have also 
demonstrated that stationary bound states of more than 
two dark solitons cannot exist. (We have, in fact, proved 
that the solutions presented in this paper are the only so- 
lutions to the reduced scalar equation (jS)) which approach 
the stable background at infinity.) 

Apart from the Ginsburg-Landau equation, Eq.Q oc- 
curred as a stationary reduction of several conservative 
nonlinear evolution equations modelling some other phys- 
ical situations. We have already mentioned the MSTB 



(Montonen-Sarker-TruUinger-Bishop) model of field the- 
ory Il9l | ; the corresponding equation of motion can 
be written as the complex ^^-equation with the broken 
C/(l) symmetry: 

l^tt-l^.. + m^^-^ + hr = 0- (59) 

(For the current status of the MSTB and related theories, 
see H^].) In addition, Eq.Q describes stationary solu- 
tions of the parametrically driven nonlinear Schrodinger 
equation: 

iA + l^x.-W^ + ^-h^* =0. (60) 

In fluid dynamics, H6()|l governs the amplitude of the oscil- 
lation of the water surface in a vertic ally vibrated chan- 
nel with large width-to-depth ratio lO, [jj. (If, con- 
versely, the channel is deep and narrow, Ea. (|60|l is still 
valid, but with the opposite sign in front of the non- 
linear term.) The same Ea. (|60|) arises as an amplitude 
equation for the upper cutoff mode in the parametrically 
driven nonlinear lattices [2^. It was also derived for the 
doubly resonant x*-^^ optical parametric oscillator in the 
limit of large second- harmonic detuning '2^. In all of 
these cases, the term hip* represents parametric pump- 
ing of some sort. Finally, Ea. (|60|l describes magnetisa- 
tion waves in a quasi-one-dimensional ferromagnet with 
a weakly anisotropic e asy plane, in a perpendicular sta- 
tionary magnetic field In the magnetic context, the 
htp* term accounts for the anisotropy of the ferromag- 
netic crystal. 

The two-soliton bound state solution we have obtained 
in this paper admits a transparent interpretation in each 
of the above physical situations. In the context of vi- 
brated water troughs and chains of coupled pendula, the 
Bloch-Neel complex describes a patch oscillating 180° out 
of phase with the rest of the chain or channel. A simi- 
lar interpretation arises in optical parametric oscillators; 
there, a bound state of two dark solitons represents a lo- 
calised fundamental field domain with a 180° phase dif- 
ference from the rest of the cavity. In the context of 
ferromagnetism, the stationary complex corresponds to 
a magnetic bubble, i.e. an "island" of one stable phase 
in the "sea" of the other one. 

The question of stability of the two-soliton bound 
state is beyond the scope of our current investigation. 
The answer will obviously depend on whether Eq.Q is 
considered as a stationary reduction of the Ginsburg- 
Landau, Klein-Gordon or nonlinear Schrodinger equation 
[i.e. Eqs.Q, (Eni or To illustrate the model- 

dependence of stability properties of one and the same 
solution, it is instructive to bring up the example of a 
free-standing Neel wall, Eq.lQ. Let, for example, h < ^. 
If the Neel wall is considered as a stationary solution 
of the Ginsburg-Landau equation ^ or of the Klein- 
Gordon equation then it is found to be unstable 
while the Bloch wall is stable H H IM El. On the 
contrary, both Bloch and Neel walls are stable 0| when 
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considered as stationary solutions of the parametrically 
driven NLS equation, Ea. l|()U|) . We are planning to re- 
turn to the issue of stability of the bound states in future 
publications. 

Finally, it is appropriate to mention that Eq.j^I) ap- 
peared in one more optical context, namely that of birc- 
fringent optical fibers. The vector nonlinear Schrodinger 
equation for pulses travelling in a birefringent fiber was 
derived by Menyuk [30j| : 



fdEi ^dEA Id^Ei 



\ dt 



dx 



2 dx^ 



\Eo 



El 



'-ElEle-^'^* 



= 0, (61a) 



. fdE2 ^dE2 
I — 



\ dt 



dx 



1^2 P I E2 



1 



ElE^e'^'^'^ = 0. (61b) 



The above assumption (5 = 0, ft. 7^ 0) can be justified 
in the case of the bright solitons, i.e. in the anomalous 
dispersion regime. In that case, Blow, Doran and Wood 
demonstrated the existence of bound pairs of bright soli- 
tons where each soliton is polarised along a different bire- 
fringence axis |3lj |. Later, these solutions were explicitly 
constructed by Tratnik and Sipe p3|- Whether the as- 
sumption 5 — Q, h ^ Q can be justified in the case of 
the dark solitons is an open question which is beyond the 
scope of our investigation. Here, we simply note that 
for time-independent real fields U and Eas. H62|l re- 
duce to our Eq.© with tp = U + iV. This fact was 
used by Christodoulides [33 who obtained the Bloch 
and Neel walls for the weakly birefringent fiber (under 
the above assumption). Our solution (|39|l is a bound 
state of the "dark" and "bright-dark" vector solitons of 
Christodoulides. 



Here S is proportional to the difference of group velocities 
of the fast and slow linearly polarised modes (whose en- 
velopes are described by Ei and E2), and h measures the 
mismatch of the corresponding propagation constants. 
Equations H61a|) - (j61bp are written in a frame moving with 
the average of the group velocities; the choice of coeffi- 
cients corresponds to the fiber in the regime of normal 
dispersion. If one assumes that the difference of the two 
group velocities is so small that it can be neglected, while 
the difference of the propagation constants is nonnegligi- 
ble (though possibly small), then the substitution 

El = C/e-^e^+i)*, ^ l/e'C'-i)*, 



takes Eq.lED to 



.dU Id^U 



dt 



2 dx^ 
1 



\U\' 



-:\V[ 



u 



.dV 



dt 2 dx^ 
1 



^V^U* + {l + h)U = 0, 



^U^V* + (1 - h)V = 0. 



(62a) 



(62b) 
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